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Abstract 

We study the pickup ion dynamics and mechanism of multiple reflection and acceleration at 
the structured quasi-perpendicular supercritical shock. The motion of the pickup ions in the 
shock is studied analytically and numerically using the test particle analysis in the model 
shock front. The analysis shows that slow pickup ions may be accelerated at the shock 
ramp to high energies. The maximum ion energy is determined by the fine structure of the 
electro-magnetic field at the shock ramp and decreases when the angle between magnetic 
field and shock normal decreases. Evolution of pickup ion distribution across the nearly- 
perpendicular shock and pickup ion spectrum is also studied by direct numerical analysis. 



1 Introduction 



Interstellar pickup ions play an important role in the physics of the outer helio- 
sphere. They may modify both the large scale characteristics of the solar wind it- 
self [1,2] and smaller scale solar wind structures such as interplanetary collisionless 
shock waves [3]. On the basis of the studies of the low-energy cosmic rays it was 
proposed that the cosmic ray anomalous component originates from the interstellar 
pickup ions, accelerated at the quasi-perpendicular shocks [4,5]. 

The observations, made by Ulysses at 4.5 AU [6], have revealed accelerated inter- 
stellar pickup ions across the forward shock in the corotating interaction regions. It 
has been found that the injection efficiency for these pickup ions exceeds that one 
for solar wind ions, and that the accelerated pickup ions have power-law energetic 
spectra in the solar wind frame [6] . Earlier investigations of the quasi-perpendicular 
cometary bow shocks have observed water-group pickup ions with speeds of sev- 
eral times the solar wind speed v u at and downstream of the shock[7,8]. 

Process of pickup ions production and features of the pickup ion distribution in the 
solar wind are studied quite well. Neutral atoms and molecules, penetrating from 
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the local interstellar medium or escaping from comets, are ionized by photoioniza- 
tion, electron impact or charge exchange with the solar wind. Under influence of 
the solar wind electro-magnetic field the ions compose ring-beam distribution in 
velocity space. Then they are scattered by ambient and excited Alfvenic fluctua- 
tions to form spherical shell distribution, centered approximately at the solar wind 
velocity with a radius of about solar wind speed V u [9,10]. 

Accelerated pickup ions were observed at quasi-perpendicular shocks, which are 
not able to accelerate incident thermal ions and pick up ions by the standard dif- 
fusive shock (Fermi) acceleration mechanism [11]. The mechanism requires ions 
which can penetrate the shock front in both directions. At quasi-perpendicular 
shocks it may occur when the ion velocity near the shock ramp is higher than a 
definite threshold speed, which significantly exceeds the solar wind speed [11,12]. 
Hence, in order that pick up ions may be accelerated by the diffusive mechanisms, 
some preacceleration mechanism at quasi-perpendicular shocks must exist. In [3], 
a detailed discussion was presented of the following idea for the ion acceleration 
at a perpendicular shock proposed by [13]: if the ion encounters the shock with 
normal kinetic energy much smaller than the electrostatic potential e(f at 

the shock, and the upstream Lorentz force is directed toward the shock, then the 
ion finds itself trapped between the shock potential and the Lorentz force. Such 
ion is multiply reflected at the shock front and during each excursion to the up- 
stream region it gains some energy until it is able to overcome the shock, when 
rriiV^/2 > eipa or the Lorentz force in the normal direction exceeds the elec- 
trostatic force in the ramp. This idea for pickup ion preacceleration at the quasi- 
perpendicular shock has been developed further by several authors. Analytical es- 
timates show that at the perpendicular shock the maximum energy gain for trans- 
mitted ions is proportional to the ratio of an solar wind ion gyroradius to the small- 
est characteristic scale of the electrostatic potential [3]. At the quasi-perpendicular 
shock maximum energy gain occurs when the ion escapes back into the upstream 
region [14]. The multiple reflection ion acceleration gives very hard power law 
spectrum [3]. The spectrum is weakly sensitive to the shock parameters and may 
extend to energy ~ 0.5MeV. Thus the injection problem associated with Fermi 
acceleration probably, may be solved by the multiple reflection ion acceleration 
mechanism. 

In the present paper we study the details of the multiple reflected ion acceleration 
mechanism in the stationary structured quasi-perpendicular supercritical shock and 
determine the shock parameters, which control the pickup ion energy gain. We 
also consider evolution of the pickup ion distribution across the strong shock. Our 
approach differs from earlier ones in that we use field models qualitatively describ- 
ing the actual structure of the stationary fields at quasi-perpendicular super-critical 
shocks which consists of the extended foot, narrow ramp, overshoot and down- 
stream region. In section 2 we consider analytically the pickup ion motion at the 
shock front and derive the multiple reflection condition as a function of the field 
parameters. In section 3 we numerically analyze the ion motion in the stationary 
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model shock front. The analysis illustrates the analytical consideration of multi- 
ple reflection process developed in section 2 and provides additional information 
about the process. The method also allows diagnostics of pickup ion distribution at 
the shock and of pickup ion spectra before the ramp and in the far upstream and 
downstream regions. 



2 Pickup ion dynamics in the shock front 



When studying the pickup ion dynamics in the shock front we adopt the usual 
assumptions that the shock is one-dimensional and stationary. In doing so we do 
not consider the effects of, for example, possible rippling of the shock surface or 
interaction with waves, assuming that the stationary electric and magnetic fields in 
the shock front determine the ion behavior (see, however, discussion in sec. 4). 

We shall work in the normal incidence frame (N), where the upstream plasma ve- 
locity is along the shock normal. Let us choose the coordinates in such a way that 
the normal is along x axis which is directed toward the downstream region (in [14] 
x axis was directed toward the upstream region), upstream and downstream mag- 
netic fields are in xz plane, and the noncoplanarity direction is along y axis. Then 
the ion motion is governed by the following equations of motion: 

m,iV x = e(E x + v y B z - v z B y ), (1) 
rriiVy = e{E y + v z B x - v x B z ), (2) 
miVz = e(v x B y - v y B x ), (3) 

where B y , B z , and E x = —d(fi/dx depend only on x, while B x = const, and the 
motional electric field E y = V U B U sin 9 = const (where V u is the upstream flow 
speed). Here subscript u denotes asymptotically homogeneous upstream param- 
eters, and 9 is the angle between the shock normal and upstream magnetic field 
vector. 

The qualitative profile of the magnetic field is known quite well. It is usually 
considered [15] to consist of the extended foot with Lf ~ 0.5(V U /Q U ) (where 
£l u = eB u /rrii is the upstream ion gyrofrequency), on which the magnetic field 
B z gradually increases by the amount of < B u . It is followed by the narrow ramp 
with the width c/oj pi > L r > c/cu pe , where the main magnetic field B z jump occurs, 
and magnetic overshoot and probably large downstream magnetic field oscillations. 
The noncoplanar magnetic field component B y is always small relative to B z , and 
negligible upstream. Substantial component of pickup ions may affect the shock 
profile and alter the typical scales (see detail discussion in [3]). Here we consider 
pickup ions as a low density test particle population. 

The electric field profile is known much worse. It is distributed over the whole 
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shock front including foot, ramp, and overshoot, and penetrates into the down- 
stream region, so that only a part of the total cross-shock potential is applied at 
the ramp.The qualitative picture of the quasi-perpendicular shock front is shown in 
Figure 1. 

An ion with a low kinetic energy in x direction (rriiV x /2 << e<j)) is unable to over- 
come the electrostatic potential at the ramp and is reflected back to the upstream 
region. If v y B z > 0, upstream Lorentz force returns it to the ramp again. In this 
way the ion becomes trapped near the ramp and quickly oscillates in x direction 
until it escapes upstream or downstream region [3,14]. 

It does not seem possible to solve (and even analyze quantitatively) the ion equa- 
tions of motion (l)-(3) in the general case of Figure 1. Instead we shall analyze 
them for the case of surfing ions, which are assumed to be trapped in the ramp 
vicinity. To do so we make an assumption (verified aposteriori) that these ions os- 
cillate quickly in x direction, while two other velocity components vary slowly on 
the oscillation period (such behavior can be expected since strong ion acceleration 
is possible only when E y acts on the ion for a substantially long time). In this case 
separating fast and slow motion, one has: 

mux = -e-j— , (4) 

ax 

rriiVy = e(E y + v z B x ), (5) 

mvz = -ev y B x , (6) 

where (4) describes fast oscillations along x, while (5) and (6) describe slow motion 
in yz plane. The effective potential 

0eff = - V y Ay - V Z A Z , (7) 

r dA y r dAz ™ 

weakly depends on time via slow time-dependence of v y and v z . Further simpli- 
fication can be achieved by consideration a model profile, where B = B u and 
E = (0, E y , 0) before the ramp, and Taylor expanding up to the first order 

By = (^)\ x=0 x, (9) 
B z = B uS in6+(^)\ x=0 x, (10) 
Ex = (^)\ x=0 x, (11) 

in the vicinity x > of the upstream edge of the ramp (x — > — oo and x — > 
+oo corresponds to the asymptotically homogeneous upstream and downstream 
regions, respectively). This approach implicitly assumes that the trapped ions do 
not penetrate the ramp deeply. In this case the effective potential takes the following 
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form 



— v y B u sin Ox, x < 0, 

\ Xx 2 — v y B u sin 0x, x > 



4 ff = , ? „ ;' (12) 



where 



dE x dB z dB y 

A= — -j + u *-r^' ( 13 ) 

ax ax ax 

where the derivatives are taken at x = 0, and eff = at the upstream edge of the 
ramp x = 0. 

Estimating typical dE x /dx ~ A<f)/L 2 ~ Q.hrriiV 2 /2eL 2 r anddB z /dx ~ AB Z /L r ~ 
B u /L r , where L r is the ramp width, one finds that the term dE x / dx in A dominates 
unless v y /V u > (V u /Q u )/L r = M A {c/ uo pi ) / L r , where Ma is Alfven Mach num- 
ber. 

Equations (4) and (12) describe oscillations in a potential well (if A > 0, which 
is typical for the quasiperpendicular shock front), with the potential minimum of 
4>mm = —{v y B u sin 0) 2 /2X at x = v y B u sin0/\. The equations of motion are easily 
solved as follows. 

In the region x < one has (the initial condition is v x = v x0 - x-component of the 
reflected ion velocity before its excursion to the upstream region, x = 0): 

v 2 x = v 2 x0 + 2v y Vt u x, (14) 

while in the region x > the solution has the following form: 

x = x m + A sin(c<jt — ip ), (15) 
v x = uA cos(u>t — ipo), (16) 

Vo = ^sin ( ^ 2 + Am ^ ye)1/2 , (IV) 
u = (e\/m t y/ 2 , (3 = v y B u , x m = 0/\, A = ^ + ^f^^ , (18) 

where v x0 is a slowly varying function of time. 

When v y and v z vary slowly, the integral § pdq = § miV x dx over a closed trajectory 
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is adiabatically invariant. Direct calculation gives 

/ = I\ + J 2 , h = f mv x dx, h = f mv x dx, (19) 

Jx<0 Jx>0 

/, = l^f , (20) 

m^o; 2 2 7T . /3 

■(/? + Am^o/e^- + arcsm 



XrriiV x0 



2(/? 2 + Am i ^ /e) 1 /2 J 



const. (21) 



It is easy to analyze the consequences in the limiting cases. If Xm^^/P 2 <C 1 
(very low initial velocities at the upstream edge of the ramp), I 2 3> h and 

n3 2 

= const, (22) 



\3/2 m y 2 



which corresponds actually to the case v x0 — > 0. In the opposite limit XrriiV 2 /(3 2 3> 
1 the upstream part dominates l x 3> I 2 and 

h ~ = const? (23) 

which immediately gives v x0 oc f^ 3 [14]. 

Maximum kinetic energy of the oscillations is determined by the energy conserva- 
tion: 

(^W = ^ + (V? u sin£) 2 /2A, (24) 
and remains relatively small with the dependence weaker than v x oc v y . 

Roughly estimating, the ion trapping ceases and it is transmitted downstream if it 
crosses the middle of the ramp. Hence the additional trapping condition is miV 2 /2 < 
e<f)(x = L r /2), or 

mVp^Vy ^ 2 /3 , VyB u L r ^ e(f) 

2 V + 2 < 2 ' ^ 

where v is the initial ion x velocity, v y0 is its initial y velocity, and is the cross- 
ramp potential. Taking into account that v is small and assuming that dE x /dx is 
the dominant contribution in A, one finds the following estimate for the maximum 
v y during the trapping: 



K) C* a ap) = _ ( ^ ) . (26) 

B u L r B u L r \B u L r v y o t 
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The second term in the last equation is a correction to the step function based esti- 
mate for (^) max [3]. 

On the other hand, the solution of the slow equations with the initial condition 

v y = v y o, v z = v z0 is 



and 



v y = v y o cos(fl u cos 6t) + (v 2 o + V u tan 6) sin(fi M cos 6t) , (27) 
v z = (v z o + K tan 9) cos(f2 u cos 9t) — v y0 sin(f2 u cos 6t) — V u tan 6, (28) 



(«v)S° = [v 2 y0 + (v z0 + V u tan 9) 2 ]^. (29) 



If (vy)^ 1 ™^ < (v y )^\ the ion is not transmitted downstream but remains trapped 
and its v y and v x decrease until it escapes upstream. 

It is difficult to determine precisely the moment when the ion escapes upstream. We 
shall estimate the escape conditions assuming that the escape itself occurs when 
\v x \ and v y come back (decrease again) to their initial values. In that case v z can be 
easily determined using energy conservation in the de Hoffman-Teller frame (HT, 
in which the plasma upstream velocity is along the upstream magnetic field). The 
transformation rule between the two frames is 

vi N) =v[ HT \ v^ = vf T \ (30) 
vi N) =v { z HT) -V u tand. (31) 

Since E x is negligible upstream, the total HT potential is zero in this region and 
energy conservation gives (v^ HT ^) 2 = const, which in our case results in the fol- 
lowing estimate of the escape velocity: 

^.escape) = _ y {N) _ ^ ^ ^ ^ 

It is worth mentioning that in this picture the ions experience reflection in which 
their noncoplanar velocity v x does not change, while others (v x and v z ) change their 
sign in HT. 

It should be mentioned that the above analysis is approximate (since the shock pro- 
file is treated in a specific model) and applies only to those ions which are trapped 
and subsequently either transmit downstream and do not appear at the ramp any- 
more or escape upstream. The analysis is not able to catch trajectories which do 
not meet these assumptions and we have to consider them numerically in the next 
section. We also do not consider here the ions which experience the shock drift 
acceleration (cf. [3]). 
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3 Numerical analysis. 



In order to illustrate the above theoretical analysis and study the features, which 
cannot be studied analytically, we perform a test particle numerical analysis of 
pickup ion trajectories in a model shock front, which is taken to resemble the ob- 
served shock profiles, as in Figure 1. The corresponding model analytical form for 
the fields is taken as proposed by [16]: 



B y (x) = A _ x „ (34) 



B x = B u cos 9 = const, (33) 

cB x dB z [x) 

AimeV u dx ' 

B z (x) = 5 U sin 0(1 + 

, 3(x + D f -3D r ), Rr-Rf^ , 
-^—(l + tanh -i X r -i) + ^-^(1 + tanh — ) 

2(x-D )\ 



+ (R - R r ) exp( 



Rd-Rr n , . , 3(s-A L -At) 

+ o (l+tanh ), (35) 

,_, , . 1 <ip era 1 <iB 2 . x 2 . 

0*0 = -f^ - n — -r^ - -7=7- ex P -^2 ( 36 ) 

ne ax 87rne ax V7ra e dj 
E y = Vu * Bu = const, E z = 0. (37) 

where d e is the scale of the additional Gaussian electric field variation, Df, D r 
and D Q control the foot, ramp and overshoot thicknesses respectively, p e >xx oc n 7e , 
n ~ B z . The suggested form of B z describes qualitatively the front structure: foot 
at —Df < x < —D r /2, ramp —D r /2 < x < D r /2, and overshoot D r /2 < x < 
(D a + Dd) (it should be noted that here for convenience x = corresponds to 
the middle of the ramp). The noncoplanar component B y is in agreement with the 
results of [17,18]. The suggested analytical profile of the cross-shock electric field 
consists of the well-known hydrodynamical part (first two terms) for polytropic 
electrons, and additional Gaussian to take into account the penetration of the field 
in deep downstream and foot and provide the observed values of the total cross- 
shock potential drop. For the specific needs of the present paper the parameters 
were chosen as follows: Rf = 1.5, R r — 5, R Q — 6, Rd = 3.2, Df = 0.3(V^/fi u ), 
D r = 0.01(V U /Q U ) (this parameter was varied, see below), D Q = 0.1(V u /il u ), 
D d = 0.15(V U /Q U ), and 7 e = 2. We chose to examine a high-Mach number Ma = 
7.5. 

The initial pickup ion distribution was assumed to be the following shell distribu- 
tion [19,20, ] 

/(v) = ^*(|v-V u |-K). (38) 
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at the upstream edge of the foot. 



Pickup ion trajectories are traced through the described shock front starting at the 
upstream edge of the foot x = —Df. Figure 2 presents the ion trajectories pro- 
jected onto the rry-plane at the quasi-perpendicular shock with 9 = 70°. A variety 
of pickup ion trajectories are seen: 1) ions which are directly transmitted down- 
stream, 2) ions which cross the ramp from downstream to upstream once or several 
times, 3) ions which are reflected in the ramp only once, 4) ions which are accel- 
erated by multiple reflection mechanism and then escape either to downstream or 
to upstream, and 5) multiply reflected ions which rotate around the ramp until they 
escape upstream. 

Next several figures provide more qualitative details about the multiple reflection 
(trapping) process. Figure 3 shows the trajectory of an ion which is transmitted 
downstream after several reflections. It is seen that at each next encounter with the 
ramp the ion penetrates into the ramp more deeply, while both v x and v y increase, 
until the ion crosses the ramp at the maximum of v y and drifts downstream having a 
high gyration velocity. The initial ion energy is about rriiV^/2, and it is accelerated 
to about four times higher downstream energy in four cycles of reflection. 

Figure 4 shows the trajectory of an ion which escapes upstream after making a 
full cycle of multiple reflections, during which v y increases and decreases again. 
Figure 4a shows that the ion trajectory does not cross the middle of the ramp. From 
Figure 4b one can see that the amplitude of oscillations of v x is roughly proportional 
to v y , that is, the adiabatic approximation works quite well. It is seen also that the 
reflection-escape process is indeed almost specular in v x — v y plane: v y is the same 
at the entry and escape points, while v x changes its sign. As can be seen from 
Figure 4c, v z in the last escape point closely corresponds to (32). 

Figure 5 shows the trajectory which is not covered by the above theoretical anal- 
ysis. The ion is trapped near the ramp for some time and afterwards it becomes 
trapped around the ramp, making several large amplitude gyrations and crossing 
the shock front back and forth. Eventually it escapes upstream (high negative final 
v z in Figure 5c) having a substantial gyration velocity. 

Figure 6a shows which part of the initial pickup ion distribution undergoes multi- 
ple reflection (for perpendicular shock geometry). In Figure 6b we show the same 
distribution of incident pickup ions at the upstream edge of the ramp, where it is 
already strongly disturbed. Almost all multiply reflected ions are taken from the 
low v x part of the distribution, in agreement with previous theoretical observations 
[3,14]. Most of them have substantial positive v y at their entry to the ramp, as is 
expected according to the analytical consideration. 

The next several figures present the evolution of the pickup ion distribution across 
the structured shock front. In Figure 7 we present the pickup ion distributions ob- 
tained numerically by tracing the initial shell distribution (38) across the perpen- 
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dicular shock. The distributions are plotted for several positions: (a) before the foot 
at x = —0.6(V U /Q U ), (b) at the upstream edge of the ramp x = —0.005(1^/^), 
(c) at the downstream edge of the ramp x = 0.005(V u /tt u ), and (d) far downstream 
at x = 5(V U /Q U ). The ion distribution before the foot (Figure 7a) consists of the 
incident pickup ion shell and ions which are reflected in the way similar to the 
reflection of the ions from the wings of the thermal ion distribution [16]. The ion 
distribution at the upstream edge of the ramp (Figure 7b) includes these reflected 
ions and much more energetic trapped accelerated ions, which are seen also just 
behind the ramp in Figure 7c (only those which crossed the ramp and are transmit- 
ted further downstream). Far downstream distribution (Figure 7d) consists of low 
energy reflected ions and high energy multiply reflected ones. Since there is not gy- 
rophase mixing in the perpendicular shock (all ions have the same downstream drift 
velocity) the downstream distribution is spatially dependent (actually periodic). 

Figure 8 present the evolution of pick up ion distribution across the nearly per- 
pendicular 9 = 80° shock. For each position both v x , v y and v y , v z projections are 
shown. The distribution before the foot (Figures 8a and 8b) differs from that one 
in the perpendicular case (Figure 7a) only by presence of ions which escaped up- 
stream (high negative v z ). The upstream edge distribution (Figures 8c and 8d) also 
shows presence of these escaping ions. The multiply reflected ions in Figure 8d lie 
on two semicircles, corresponding to the slow rotation in v y , v z plane (see (5) and 
(6)), separately for positive and negative initial v z . The most remarkable difference 
from the perpendicular case can be seen in the downstream distribution (Figures 8g 
and 8h), which shows strong phase mixing. The downstream accelerated pickup 
ions are situated on two (almost) hemispheres, corresponding to the rotation of the 
semi-circles in Figure 8f. 

The numerical analysis allows also to obtain the differential energy spectra of ac- 
celerated ions. Figure 9 shows the far upstream (Figure 9a) and far downstream 
(Figure 9b) ion energy spectra dN / de, 



where e is the dimensionless ion energy, e = (rriiV 2 /2)/ (m>iV 2 /2). The far up- 
stream spectrum consists of the dense population of incident ions (dN/de = const) 
and low density high energy component of accelerated ions which escaped up- 
stream for two different angles between the shock normal: 9 = 80° (dotted line) 
and 9 = 70° (dashed line). With the decrease of the angle the number of escape 
ions (Figure 9a) increases but their energy decreases. One can see also the drastic 
drop of acceleration efficiency with the increase of obliquity in the downstream dis- 
tribution of pickup ions, where the highest energy of accelerated ions drops from 
10 2 (in rriiV 2 /2) to slightly higher than 10 1 when the angle decreases from 80° to 
70°. 

Figure 10 shows the accelerated pickup ion distribution just before the ramp and far 




(39) 
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downstream for two values of the total cross-shock potential. As could be expected, 
the acceleration efficiency decreases with the decrease of the potential, although 
this decrease is almost not noticeable in the far downstream distribution because 
of the logarithmic scale. It is clearly seen that the distribution is nearly exponential 
dN/de oc exp(— cte) at the upstream edge of the ramp, with a 0.020 for the po- 
tential up = 0.7(miV^/2), and a « 0.033 for up = 0.5(m;\/ 2 /2). The downstream 
distributions reveal power spectra dN/de oc e _/3 , where (3 ~ 3/2 and is almost in- 
dependent of the cross-shock potential. The result is in conformity with the power 
low tail in energy produced by multiply reflected ion acceleration at strong perpen- 
dicular shock found in [3]. 

Dependence of the accelerated ion spectra on the ramp width is shown in Fig- 
ure 11, where the upstream-ramp-edge and far downstream spectra are presented 
for the ramp width D r = 0.01(KMJ (solid line) and D r = 0.03(V U /Q U ) (dotted 
line). As is expected, the acceleration efficiency drops drastically when the sock 
becomes wider. It could be expected, however, that any substructure in the ramp 
would enhance the acceleration, as can be seen from Figure 12, where the ion spec- 
tra are compared for the case of D r = 0.03(V u /tt u ) without (dotted line) and with 
a substructure (dashed line). Presence of such internal substructure plays the role 
of a narrow ramp. 

It is worth to mention that our approach allows also to study the behavior of heavy 
ions. Figure 13 presents the results of such analysis, comparing the downstream 
spectra of ions with the mass m 8 = Am p (which would correspond to He + ions) and 
Tfii = 15m p (0 + ). The initial distribution is the same pickup ion shell as for the case 
rrii = m p , considered throughout the paper. It can be seen that the acceleration is 
less efficient for heavy ions than for protons (Figures 9b and 10b). Detailed analysis 
of the heavy ion behavior is beyond the scope of the present paper and will be 
presented elsewhere. 



4 Conclusions 

We have considered the surfing mechanism of the pickup ion energization at strong 
high-Mach number quasi-perpendicular shocks. The distribution of the electro- 
magnetic field at the shock front determines the ion motion across the shock,and 
the ion reflection and acceleration processes depend on the details of the fine struc- 
ture of the shock front. We had to make a choice of the shock structure in order to 
determine quantitatively the features of the accelerated ions and their dependence 
on the shock parameters. It is clear that a deviation of the actual shock structure 
from the model adopted here would give somewhat different quantitative results, 
which cannot be predicted unless we know profiles of the fields at the shock. How- 
ever, some tendencies can be predicted on the basis of the analytical and numerical 
investigations represented in sections 2 and 3. 
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The conditions for the multiple ion reflection and the maximum energy gain are 
determined by the slope of the electric and magnetic field profiles at the ramp and 
therefore very sensitive to the ramp width. For the chosen model the scales of the 
magnetic and electric field variations are the same and ss D r . The analytical con- 
sideration gives the following estimate for the maximum downstream ion energy in 
the nearly perpendicular shock (see also [14,3]): 

where (j) — eip/ (m* , R = B zr / B u (where B zr is z component of the magnetic 
field at the downstream edge of the ramp), and D r is measured in V u /Q u . For the 
initial ion energy of S = IkeV, ramp width of D r 0.01(V u /Q u ), cross-shock 
potential <fi = 0.5, and magnetic compression ratio R = 2.5, one finds £ max ss 
lMeV, which is probably sufficient for injection into diffusive acceleration regime. 

Since the efficiency is oc D~ 2 and decreases with the increase of obliquity, the pick 
up ion injection for Fermi acceleration mechanism seems to be effective only at 
nearly-perpendicular shocks waves with very narrow ramp or when some substruc- 
ture is present which reduced the effective width of the ramp. 

The dependence of the acceleration features on the Mach number is determined 
mostly by the dependence of the shock width and magnetic compression ratio on 
the Mach number. In the absence of a satisfactory theory, which could provide these 
dependencies, the only conclusion is that the mechanism efficiency should rapidly 
decrease with the decrease of the Mach number because of the shock widening. 

The obtained model downstream accelerated ion spectra oc e~ 3 / 2 (see also [3]), 
which corresponds to f(v) oc v ~ 2 (and is in agreement with earlier studies), is 
harder than v ~ 4 that is observed. This discrepancy may be probably attributed to 
the mentioned deviation of the actual shock structure from the adopted model. 

It should be mentioned that, due to the high sensitivity of the trapping-detrapping 
mechanisms to the details of the magnetic and electric fields in the ramp, it could 
be in principle affected by even relatively weak deviations from stationarity (for ex- 
ample, presence of large amplitude waves in the ramp) or one-dimensionality (like 
rippling of the shock surface), as well as presence of any small-scale substructure 
inside the ramp. Observation evidence is not unambiguous, and the effect of such 
deviations from the model is not clear apriori and requires special study which is 
beyond the scope of the present paper. 

In summary, the analysis allowed to shed additional light on the importance of 
the small-scale structure of the shock front (in particular, high gradients of the 
electric field) for the ion acceleration processes. The test particle numerical analysis 
allowed to determine the details of ion behavior and the spatial and energetical 
distribution of accelerated ions in the shock front. 
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